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Abstract 

We investigate the dynamics of the chameleonic Generahsed Brans-Dicke model in flat FRW 
cosmology. In a new approach, a framework to study stability and attractor solutions in the phase 
space is developed for the model by simultaneously best fitting the stability and model parameters 
with the observational data. The results show that for an accelerating universe the phantom 
crossing does not occur in the past and near future. 

PACS numbers: 04.20.Cv; 04.50.-h; 04.60.Ds; 98.80. Qc 

Keywords: chameleon; Brans-Dicke; stability; attractor; distance modulus; method 



*Electronic address: 
^Electronic address: 



hosseinf @ guilan . ac . ir 



a. salehi® guilan . ac . ir 



1. INTRODUCTION 



Recently, the observations of high redshift Type la supernovea (SNe la), the surveys of 
clusters of galaxies PP-|1], Sloan digital sky survey (SDSS) [5] and Chandra X-ray ob- 
servatory |6] reveal the universe accelerating expansion. Also the observations of Cosmic 
Microwave Background (CMB) anisotropics ^ indicate that the universe is flat and the total 
energy density is very close to the critical one [8] . The observations determines basic cosmo- 
logical parameters with high precisions and strongly indicates that the universe presently is 
dominated by a smoothly distributed and slowly varying dark energy (DE) component. A 
dynamical equation of state ( EoS) parameter that is connected directly to the evolution of 
the energy density in the universe and indirectly to the expansion of the Universe can be 
regarded as a suitable parameter to explain the acceleration and the origin of DE |9]-|10]. 
Among all the cosmological models, the scalar-tensor theories have been widely used to ex- 
plain the late time acceleration of the universe and its relation to the dark energy. For a 
review see 

The authors in [33] in a separate work have studied the dynamics of the universe in 
Chameleonic Brans-Dicke (CBD) cosmological model. They perform stability analysis and 
investigate the attractor solutions of the CBD cosmology by utilizing the 2-dimensional 
phase space of the theory. In there, for the matter in the universe assumed to be perfect 
fluid with the EoS parameters, 7 = or 1/3, the analysis gives the corresponding conditions 
for tracking attractor and determines the universe behavior in the past and future. The 
authors then study the cosmological parameters such as effective EoS parameter for the 
model in terms of the stability variables. Depending on the stability parameters, the model 
predicts phantom behavior, despite the absence of phantom energy. It has been shown that 
by fitting the observational data to the model for the distance modulus, the scenarios with 
phantom crossing better fits the data. 

In an extension, in this work we study the Chameleonic Generalized Brans-Dicke (CGBD) 
cosmological model by letting the BD parameter to be function of the BD scalar field. A 
dynamical BD parameter can be favored by realizing that general relativity is recovered from 
BD model in the limit w — > 00 and large values of the BD scalar field. Also, in a different 
approach, to perform stability analysis, we simultaneously best fit the model parameters 
with the observational data using method. This enables us to find the best fitted model 
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parameters for tlie analysis of the critical points and also verification of the model with the 
experiment. 

The paper is organized as follows: Sec. two is devoted to a detailed formulation of the 
CGBD cosmological model. In Sec. three, we perform stability analysis of the model. We 
first obtain the autonomous differential equations of the model in terms of dimensionless 
dynamical variables. We then solve the equations by simultaneously best fitting the model 
parameters and initial conditions with the data and study the attracting behavior of the 
critical points in the phase plane. In Sec. four, we examine the behavior of the EoS parameter 
of the model. In Sec. five, we present a discussion. 

2. THE MODEL 

We consider the CGBD gravity in the presence of matter with the action given by. 



where R is the Ricci scalar, is the Chameleon and/or BD scalar field with the potential 
V{(f)), and uj{(j)) is BD parameter. We assume that /(0) = f^cj)'^ and V{(j)) = Vocj)^ where k, 
and 6 are dimensionless constants characterizing the slope of potential V{(j)) and /(0). The 
cosmological models with such exponential functions have been known lead to interesting 
physics in a variety of context, ranging from existence of accelerated expansions [35j to 
cosmological scaling solutions [36]-[33. In particular, the exponential forms of /(0) and 
V{(f)) are motivated by chameleon models [32] and also from stability considerations 
In addition, attractor solutions with exponential functions may lead to cosmic acceleration 
for natural values of model parameters [39]. Unlike the usual Einstein-Hilbert action, the 
matter lagrangian Lm is modified as /(0)Lm, where /(0) is an analytic function of the scalar 
field. The last term in the action brings about the nonminimal interaction between matter 
and the scalar field. Variation of the action (IT]) with respect to the metric tensor g^^ gives. 




(1) 
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2i/ + 3//-^ = -^(^-2i/^-^^-^-'^. (4) 

(p (j) 2 (j)^ (j) (j) 

On the other hand, the variation of the action ([T]) with respect to the spatially homogeneous 
scalar fields 0(t) gives, 

^ ^ 3 + 2u;(0) 3 + 2w(0) 2(3 + 2w(0)) 3 + 2w(0)' ^ ' 

where prime " ' " means derivative with respect to 0. We assume that the universe is filled 
with the barotropic fiuid with the EoS, pm = IPm- From equations ([3]), Q and (|5]), one can 
easily arrive at the modified conservation equation, 

(p„/) + 3^(1 + 7)Pn^/ = ^(1 - 37)Pm/. (6) 

In the next section we study the stabihty analysis of the model in the phase space. 



3. PERTURBATION AND STABILITY ANALYSIS 



In this section, we study the structure of the dynamical system via phase plane analysis, 
by introducing the following dimensionless variables. 



X^ = ^ C = — ^7^ = — (7) 



Using equations ([3])-(j5]), the evolution equations of these variables become. 



-3(1 + 7) , r/€(l- 37)- 20 /3a;oe°^ , uoe^"" 

6 + 4a;oe"^^ 



x' = x{ ^ ' " + r - . : j c + (s) 



3(l-37)(2-/.) 3(1 + 7) . 2 , 3^ - 6 ^ 
U(3 + 2u;oe°^) 2 3 + 2u;oe°^^ ^ 



aa;oe°^ c^oe"^ 2 . 3(1 - 37)(2 - /t) , 3(1 + 7)^ 
6 + 4w^ 2 ^ ^ 4(3 + 2u;oe^ 



c = c{-3 - 3c - ^^^c + + r:^. + -^u^ o) 

3 + 2woe"^ '' ^ ' 2(3 + 2u;oe"^) ' 3 + 2woe"^^ ' 



3(5-6 21 , 3(l-37)(2-/t) 2 , 6^-12 2 



, = _ ^W!^ 2 + r 3(l-37)(2-.) 3(1 + 7) 2 . 

' 2 ^ 6 + 4u;oe"^^ 2 ^ ^ 2(3 + 2u;oe°^) 2 ^ ^ 

35 — 6 9, 

}, 



3 + 2woe"^ 

where prime " ' " here and from now on is taken to be derivative with respect to = ln{a). 
The parameter a in the above equation is defined by a = ^ which gives u = Uqc""^ . Also, 
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the Friedmann constraint equation (|3| in terms of the new dynamical variables becomes 



X -C + — -r/ = 1. (11) 

In term of the these variables we also obtain, 

«^oe"^ ^oe"^ 2 . 3(l-37)(2-/.) 3(1+7) . 3^-6 , 

i72 ^ ^6 + 4u;oe"^^^ 2 ^ U(3 + 2a;oe°^) 2 3 + 2u;oe"^^^ ' 



Using the Friedmann constraint (11), the three equations (|8j)-(10) reduce to the following 
two equations: 



-3(1 + 7) , ^k(1 - 37) - 20 awoe"^ ^ , Wqc"^ 



x' = + r ' - . j c + (la) 



. 3(1-37)(2-/.) 3(1 + 7) , 3^-6 . ^ + 2 _ . ^ 
^ 4(3 + 2woe°^) 2 ^-^ 3 + 2^06"^ 6 ^ 



^ ^1 '^^ 6 + 4a;oe"^^ 2 ^ ^ U(3 + 2a;oe"^^ + ^"'^^ 



+ ^^^(-l + v^-C + ^ai + ^ii^M^.2 
3 + 2a;oe°^^ ^ 6 ^ 4(3 + 2cJoe"^) ^ 



It is more convenient to investigate the properties of the dynamical system equations(13) 



and (14) than equations ([8j)-(10). We, now, obtain the critical points and study their 
stability of these points. Critical points are always exact constant solutions in the context 
of autonomous dynamical systems. These points are often the extreme points of the orbits 
and therefore describe the asymptotic behavior of the system. In the following, we find 
fixed points by simultaneously solving x' = and C'- Substituting linear perturbations 
x' x' + ^x' ^ C C + about the critical points into the two independent equations 



(13) and (14), to the first orders in the perturbations, one yields two eigenvalues Aj(i = 1, 2). 
Stability requires the real part of the eigenvalues to be negative. 

In the following, the above equations can be solved to give fixed points (critical points) 
for two scenarios of the matter field, 7 = and 7 = 1/3. Both critical points and eigenvalues 
in our model are highly nonlinear and depend on the stability parameters k, 5 and model 
parameters a and uq. In addition, the expressions for the critical points and eigenvalues are 
long and cumbersome such that it is not easy to find out under what conditions the critical 
points are stable or unstable. 
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In a different approach we solve the above equations by best fitting the stabihty and model 
parameters and initial conditions with the observational data for distance modulus using the 

method. This helps us to find the solutions for the above equations and conditions for 
the stability of the critical points that are physically more meaningful and observationally 
more favored. Next, we examine this idea of simultaneously fitting the model and solving 
the field equations. 



3.1. Best fitting tiie stability parameters and initial conditions 

The difference between the absolute and apparent luminosity of a distance object is 
given by, fi{z) = 25 + 5 log^Q ^2,(2;) where d^i^z) is the Luminosity distance quantity. If we 
consider a possible variation of the effective gravitational constant, G^ff, in the Universe 
acceleration rate history, an additional term made of with the ratio between the value 
of effective gravitational constant at any redshift and the same quantity evaluated at the 
present has to be contributed, |10] 

^{z; /3, L; A) = 25 + 5 log^o d,{z) + ^ log "^'"^''^'^'^ (15) 

where 

G^ffiz; P, L; A) =Gm{1 + ^/^ ( /| p ) ^ 

with Gn is the Newtons gravitational constant. The intrinsic model parameters, i.e., the 
coupling constant P and the interaction length L oc m^^ can be constrained with the fitting 
procedure. The length A oc (wavelength k) for the SNe la redshift range is of order 
10~^. As explained in |10], a time- varying gravitational constant can affect light curves from 
SNe la. However, the correction coming from it is negligible when compared to the usual 
distance modulus expression. 

From numerical computation one can obtain H{z) which can be used to evaluate fj^{z). 
To best fit the model for the stability and model parameters S, k, a and Uq and the initial 
conditions x(0), C(0), H{0) with the most recent observational data, the SNe la, we employe 
the method. In addition, from [H] we use the best fitted values obtained for the parame- 



ters (3 and L in calculating the last term of the equation (16). We constrain the parameters 
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and also the initial conditions by minimizing the function given by 

r fhef^ \x ,^ ^, , , ■.,fn\ /"^'n^ UtnW ,,obsl2 



XsNei^^ K, a, Wo, X(0), C(0), H{0)} = > 2 ' (16) 

— (Jf 



1=1 

.the „„J ,,obs 



where the sum is over the SNe la sample. In the relation (16), /i* and /i° are the distance 
modulus parameters obtained from our model and from observation, respectively, and a is 
the estimated error of the ^f"^. From numerical computation in table I, we show the best 
fitted model parameters for 7 = and 7 = 1/3. 



TABLE I: Best fitted values for the Stability parameters, model parameters and Initial conditions. 



parameters 


5 K 


a 


UJo 


x(o) c(o) ^(0) 


• 


7 = 


-3.1 22 


-1.3 


9.7 


-0.3 0.2 0.91 


543.1628893 


7 = 1/3 


-0.04 - 


-1.78 


6.2 


-0.3 0.1 0.91 


552.4495936 



Fig. 1) shows the constraints on the parameters S, k, a and and the initial conditions 
X(0), C(0) and H{0) at the 68.3%, 95.4% and 99.7% confidence levels in both cases of 
7 = 0, 1/3. Note that in case of 7 = 1/3, since in equation (12) the third term vanishes, no 



constraint on k is expected, as can be seen in Fig.l. 




CO 




-2.3 -1.3 -0.3 



Fig. 1: The constraints at the 68.3%, 95.4% and 99.7% confidence levels 

from Sne la for the parameters 6, n, a and Uq 

and the LC.s x(0), C(0) and H{0) in case of: left & middle) 7 = right) 7 = 1/3 
Alternatively, we can plot the likelihood for the model parameters in both cases 7 = and 

7 = 1/3 ( Figs. 2 and 3). 
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Fig. 2: The 1-dim and 2-dim likelihood distribution for parameters k, 6, Uq and a in 7 = case. 





Fig. 3: The 1-dim and 2-dim likelihood distribution for parameters Uq and a in 7 = 1/3 case. 



The distance modulus, fi{z), plotted in Fig. 4, in both cases 7 = 0, 1/3 are fitted with the 
SNe la observational data for the model parameters and initial conditions using method. 
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0.5 1 1.5 0.5 1 l.S 



z z 

Fig. 4: The distance modulus fi{z) plotted as function of redshift for the best fitted parameters 

in both cases 7 = 0, 1/3. 

In the following we investigate the stability of the model with respect to the best fitted 
model parameters and for the two specific choices of the EoS parameter for the matter in 
the universe, i.e. 7 = and 7 = 1/3. 

3.2. Stability of the best fitted critical points and phase space 

Solving the stability equations for the best fitted model parameters we find fixed points 
with their stability properties as illustrated in tables II and III, for 7 = and 7 = 1/3 
respectively. 



TABLE II: Best fitted critical points for 7 = 



points (x, C) 


Stability 


PI (0, 0) 


unstable 


P2 (0.78,0.37) 


stable 



P3 (-0.78, .37) stable 



From the above tables we see that, by best fitting the stability parameters, the critical 
points P2 and P3 for 7 = and P4 for 7 = 1/3 are stable and the others are unstable. 
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TABLE III: Best fitted critical points for 7 = | 



points 


ix, C) Stability 


PI 


(0, 0) unstable 


P2 


(1, 0) saddle point 


P3 


(—1,0) saddle point 


P4 


(0, 4.3) stable 



In Fig. 5)left), for 7 = 0, the trajectories leaving the unstable critical point PI in the past 
in the phase plane is shown going towards the stable critical point P2 or p3 in the future. 
The best fitted model parameter trajectories is also shown by red color dashed trajectory). 

In Fig. 5)right), for 7 = 1/3, there are three unstable critical points including two 
saddle points, PI and P2. With the given initial conditions for x(0) between -1 and 1, the 
trajectories (blue and red curves) leaving the unstable critical point PI in the past in the 
phase plane is shown moving towards the stable critical point P4 in the future. For the 
initial conditions less than -1 or greater than 1, then the trajectory shown by green curve 
leaving the unstable critical point PI moving towards the saddle points P2 and P3 and 
leaving these points towards infinity. 




X X 

Fig. 5: The attractor property of the dynamical system in the phase plane for 

left) 7 = and right) 7 = 1/3. 
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4. COSMOLOGICAL PARAMETERS: EOS PARAMETER 



In order to understand the behavior of the universe and its dynamics we need to study 
the cosmological parameters such as EoS parameter. We have aheady fitted our model with 
the current observational data for the distance modulus. The EoS parameters analytically 
and/or numerically have been investigated by many authors for variety of cosmological 
models. Applying stability analysis and simultaneously best fitting the model with the 
observational data using method gives us a better understanding of the universe state. 
The effective EoS parameter is defined by Wejf = —1 — where is given in terms of 



best fitted new dynamical variables in equation (12). 

The effective EoS parameters obtained in both cases, 7 = 0, 1/3, do not exhibit phantom 
crossing behavior for different stability parameters including the best fitted ones. In case 
of 7 = 0, the EoS parameters for different model parameters including the best fitted 
one corresponding to the unstable critical points PI in the past start from Weff = 1, and 
approaches w^ff = —0.5 in the future corresponding to the stable critical point P2 or P3. 
The current effective EoS parameter for different stability parameters changes and the best 
fitted value is about —0.6. Numerical solutions also reveal that the effective EoS parameter 
begins to oscillate around Weff = —0.5 at late time, as shown in the left panel of Fig. 6. 
The amplitude of the oscillations are decreasing with respect to ln{a). 

In case of 7 = 1/3, the effective EoS parameters for different model parameters including 
the best fitted one corresponding to the unstable critical point PI start from Weff = 1 in the 
past, and approach Wejf = 0.48 in the future corresponding to the stable critical point P4. 
The current EoS parameter for different stability parameters changes and the best fitted 
value is —0.79. One of the trajectories (green curve) starts from PI passes P2 or P3 and 
finally tends to infinity in the future. 
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-10 -5 5 10 -10 -5 5 10 

-ln(l-i-z) -ln(l+z) 

Fig. 6: The effective EoS parameter Wefj, plotted as functions of redsliift 
for left) 7 = and right 7 = 1/3. 



5. DISCUSSION 

This paper is designed to study dynamics of the CGBD cosmology by using the stability 
analysis and the 2-dimensional phase space of the theory. The potential V{(j)) and coupling 
scalar function /(0) in the model is assumed to be in power law forms in the phase space. 
The matter lagrangian in the model is regarded as a perfect fluid with two kinds of EoS 
parameters, i.e. 7 = 0,1/3. In a different approach in stability analysis, here, we solve 
the system of autonomous differential equations by best fitting the model parameters and 
also the initial conditions with the observational data for distance modulus. Therefore all 
the critical points with the stability conditions are presented in the model are physically 
reliable and observationally verified. By stability analysis, the critical points in the model 
for the two kind of matter fields are evaluated and shown in Figs. 5. In comparison with 
our previous work in the stability analysis of CBD model where we assume a constant 
BD parameter, and no best fitting of the stability parameters is performed, in the current 
work, BD parameter is set to be an exponential function of the scalar field and best fitting 
of the model and stability parameters is also implemented. 

We then study the cosmological parameters such as effective EoS parameter, Weff for 
the model in terms of new dimensionless dynamical variables. The model predict a late 
time universe acceleration expansion. While in CBD model phantom behavior occurs for 
the universe, in this work the best fitted effective EoS parameter does not display crossing 
the divide line, in both 7 = 0, 1/3 cases. It can be verified that the CGBD model for the 
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constant BD parameter {a = 0) become CBD model. In comparison of the observational data 
for the EoS parameter of the matter, 7 = 0, 1/3, the results show that the current values for 
Weff — —0.7 corresponding to the best fitted model parameters confirm the new findings. 

The oscillating behavior of the effective EoS parameter in case of 7 = 0, is due the 
spiral property of its attractor. A phase coherent or spiral attractor is one of the typical 
nonhyperbolic chaotic regimes that decays according to an exponential law with a decrement 
which is defined by the phase diffusion coefficient |12|. It has been shown numerically and 
experimentally that spectral and correlation properties of these attractors can be adequately 
described by the model of a random process of the harmonic noise type |13]. The damped 
oscillation for the effective EoS parameter corresponding to a universe undergoing late-time 
oscillating acceleration until it reaches its stable state in the future. 
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